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Abstract. Understanding the statistical properties of the aperiodic planar Lorentz gas stands 
as a grand challenge in the theory of dynamical systems. Here we study a greatly simplified but 
related model, proposed by Arvind Ayyer and popularized by Joel Lcbowitz, in which a scatterer 
configuration on the torus is randomly updated between collisions. Taking advantage of recent 
progress in the theory of time-dependent billiards on the one hand and in probability theory on 
the other, we prove a vector-valued almost sure invariance principle for the model. Notably, the 
configuration sequence can be weakly dependent and non-stationary. We provide an expression 
for the covariance matrix, which in the non-stationary case differs from the traditional one. We 
also obtain a new invariance principle for Sinai billiards (the case of fixed scatterers) with time- 
dependent observables, and improve the accuracy and generality of existing results. 
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I. Introduction 

Recall that having infinitely many fixed scatterers in a periodic configuration in M. 2 and a particle 
moving in the exterior of the scatterers, the particle being in free motion up to elastic collisions 
with the scatterers, is a model often called the planar periodic Lorentz gas. For example, if the 
configuration is invariant under the translations of Z 2 , the model corresponds to Sinai billiards on 
the two-dimensional torus. Under certain assumptions, including that the free path of the particle 
is uniformly bounded from above and below, many statistical limit results are accordingly known 
to hold true. A planar aperiodic Lorentz gas is obtained by relaxing the periodicity assumption 
on the scatterer configuration. For example, one can begin with a periodic configuration, shift 
each scatterer by a small amount, independently of the others, and fix the resulting configuration 
for good. Understanding the statistical properties of the particle trajectories in this case is an 
outstanding problem in modern dynamics. A source of major difficulties is the phenomenon or 
re-collisions: The billiard particle could hit the same scatterer infinitely many times and, because 
of that, the scatterer configuration seen by the particle at a given time depends in a complicated 
way on the history of the billiard trajectory. On the other hand, if the configuration has temporal 
randomness — say it is refreshed randomly and independently of the past after each collision — the 
situation is simpler. Yet, little if anything is presently known even in that setting. The motivation 
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of this paper is to improve the state of the affairs. To that end, we are going to analyze a model 
of Sinai billiards with random scatterers, in which the configuration is randomly updated between 
collisions. 

Consider, then, the following random billiard table on a two-dimensional torus obtained by 
identifying the opposite sides of the unit square. First, a disk of radius R is placed on the torus 
with its center located at the corner of the disk, after which it is fixed for good. We will refer to 
this disk as the "gray disk" in the future. Next, another disk of radius R, called the "white disk", 
is placed on the torus at a random location so that its center is within distance e from the center 
of the square; see Figure 1. We will assume throughout that the "no-overlap" conditions 



hold. Given such a triple (R, R, e), the strict inequalities guarantee a uniform positive lower bound 
on the distance between the disks, regardless of the centering of the white disk. 

The rules of the dynamics on the table are as follows. Consider a particle traveling with unit 
speed in the complement of the disks on the surface of the torus. If the particle hits the white disk, 
it bounces elastically off the inert boundary of the disk and continues its motion with unit speed, 
as illustrated in Figure 1. If, on the other hand, the particle hits the gray disk instead, it behaves 
in exactly the same way, but the location of the white disk is refreshed randomly, so that its center 
remains within distance e from the center of the square. Finally, corresponding to the sides of the 
original unit square, there are four "transparent walls" (the opposite pairs of which are identified) 
terminating on the gray disk. If the particle meets a transparent wall, it simply passes through. 
However, if it makes a "clean pass" in that it avoids hitting the gray walls at the two ends of the 
transparent wall in question, see Figure 2 and its caption, the location of the white disk is again 
refreshed randomly. In brief, the dynamics is that of ordinary billiards, except that the location 
of the white disk is refreshed randomly each time the particle hits the gray disk or makes a clean 
pass through a transparent wall. 



Figure 1. The random billiard table. A two-dimensional torus is obtained by 
identifying the opposite sides of the unit square shown. The gray disk of radius R 
is fixed for good with its center at the corner. The white disk has a fixed radius R 
and a random center within distance e from the center of the square. The random 
centering is refreshed after each collision with the gray disk or a clean pass through 
a transparent wall; see text for details. 



max(i?, R + e) < \ and R + R + e < 



i 




Notice that the no-overlap conditions with their strict inequalities yield a uniform positive lower 
bound on length of the free path of the particle. It will be necessary to bound the length of the 
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Figure 2. A velocity vector corresponding to a "clean pass" through a transparent 
wall. A pass is clean if the segment of unit length parallel to the velocity vector and 
centered at the point of crossing does not intersect the gray disk. 



free path (between successive collisions with a disk) from above. To that end, we introduce the 
so-called "finite horizon" conditions on the geometry of the disks that we will assume to be in 
effect at all times. With the aid of Figure 3, it is easy to verify that the condition 

guarantees there is no diagonal passage along which the particle could escape without ever colliding 
with a disk. On the other hand, the condition 

R+R-e > \ 

guarantees there is neither vertical nor horizontal passage; the white disk will then, in any allowed 
position, contain the smaller centered disk shown in Figure 3. In fact, under these conditions the 
free path is uniformly bounded above. 




Figure 3. The finite horizon assumption. If the gray disk has radius R > 

the diagonal passage is blocked. Given R, the white disk blocks the horizontal and 
vertical passages if it contains the centered dashed disk of radius \ — R. 



Of course, if the above conditions are satisfied for some value of e, given R and R, they remain 
satisfied for any smaller value of e as well. The idea is that R and R are considered fixed once 
and for all and, for all sufficiently small values of e, the statistical properties of the motion of the 
particle can be analyzed. 

To the author's knowledge, the problem studied here was initially proposed by Arvind Ayyer, 
after which it got rather widely advocated inside the mathematical physics community by Joel 
Lebowitz. It was also put forward in [1], and the recent [17] mentions a reminiscent model due to 
Lebowitz. So far, no results on it seem to have been published in the literature. This is somewhat 
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surprising, given its relation to the aperiodic Lorentz gas and the fact that the stochastic stability 
properties of a similar model were studied already in [28]. For a toy version of the model, some 
theorems were obtained in [1] and [2]. See also [14,23] for related models. Below we prove a 
fairly sophisticated limit result for the model, namely the almost sure invariance principle for 
vector-valued observables. Instead of providing the precise statement here, let us formulate a 
"pre-theorem" , which captures the essence: 

Pre-theorem. There exists a Brownian motion B such that the unfolded planar trajectory 1 of 
the particle, for a typical initial condition, for a typical sequence of scatterer configurations, is 
approximated by a sample path of B. 

As the aim of the paper is to demonstrate a phenomenon, above all, we have not shied away 
from making convenience assumptions on the geometry at the expense of generality. See, however, 
Remark 3 on some immediate generalizations of no extra cost. 

The analysis below relies on several standard constructions in the theory of dispersing billiards. 
These include homogeneous local stable and unstable manifolds, together with a good understand- 
ing of their properties, among others. In particular, the random-scatterer model we consider here 
differs from the classical Sinai billiards in which the configuration of scatterers remains forever 
fixed. The excellent book [6] contains a detailed account of the classical theory and the paper [26] 
the necessary uniform time-dependent generalizations. (Also [18, 25] might interest the reader, 
although the setup there is simpler than in [26].) It is also necessary to make the observation that, 
although a billiard map is traditionally defined as the first-return map of the billiard flow to the 
solid boundary of the domain in question, the theory extends without difficulty to cover (flat or 
curved) transparent walls as cross sections [11, 12]. For the sake of clarity, we will not dwell on 
these quite involved but well-understood constructions here, instead referring the interested reader 
to the references cited. The mandatory deviation from this occurs when we introduce the billiard 
maps: They need to be defined with due care, for the presence of flat (here transparent) walls does 
have the tendency to render the billiard maps non-uniformly hyperbolic. In our case the system 
remains uniformly hyperbolic by a careful choice of the cross section. Other than that, the focus 
in the present paper is solely on establishing the estimates required to prove the desired invariance 
principle. 

An almost sure invariance principle for Sinai billiards (fixed scatterers) was first obtained for 
scalar- valued observables in [5,6,15] and for vector- valued ones in [16]. The approach in the 
present paper — which in particular applies to the setting of fixed configurations of those works 
- is based on verifying the conditions of [8] on the characteristic function of the vector-valued 
process. As the title suggests, the latter paper has been tailored for dynamical systems specifically 
with the transfer operator formulation in mind. For billiard models the spectral method [7] is 
quite cumbersome and abstract due to the well-known issue of singularities, on top of the presence 
of a contracting direction. Already the contracting direction would typically involve studying the 
transfer operator on a Banach space of distributions, as in the case of Anosov systems. Here we 
take an alternative route, providing sufficient information on the characteristic function via strong 
bounds on correlation functions. In this sense, our work is related to [24] and [5]. The fact that 
we succeed in doing so also underlines the possibility of implementing [8] in situations where the 
operator setup may be impractical or out of reach. 



Between two successive times when the scatterer configuration is updated, the particle experiences a displace- 
ment which has a natural representation as a plane vector. Summing these yields a trajectory in the plane. 
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Let us point out that an almost sure invariance principle for vector-valued observables is a sig- 
nificant improvement compared to the case of scalar- valued observables. Although it is a standard 
and straightforward trick to deduce a multi-dimensional central limit theorem, where convergence 
takes place in distribution, from the corresponding one-dimensional result using the Cramer- Wold 
theorem, the same cannot be said about almost sure invariance principles. 

How the paper is organized. In Section 2 we give a mathematical description of the model. 
Once the necessary concepts have been introduced, we present the main results of the paper in 
Section 3. There we also discuss the approach used in its proof. Section 4 is devoted to the 
technical preliminaries that are required for the analysis of the problem. In Section 5 we prove the 
main results by establishing sufficient bounds on certain correlation functions and by analyzing 
the structure of the covariance matrix. 



In this section we describe the model in full detail so that the main theorem can be formulated 
in the next section. 

2.1. The free zone condition. In addition to the no-overlap and finite horizon conditions, we 
assume the "free zone" condition 



on (R, R, e). By elementary geometrical computations, the free zone condition can equivalently be 
stated as follows: 

Given the triple (R,R,e), the white disk does not intersect any chord originating 
from a transparent wall and terminating on a piece of the gray disk at the end of 
the same transparent wall, in any possible location of the white disk. 

This is clarified in Figure 4. The role of the free zone condition is one of convenience, and it could 
be done away with. As our objective is not to be as general as possible (see Introduction), we 
choose to keep it. Under the free zone condition, it will be obvious to the reader that the billiard 
trajectory will never hit the white disk more than once between successive returns to the cross 
section M. to be defined below. 

In what follows, we will assume that the triple (R, R, e) has been fixed so that all the conditions 
introduced so far hold. 

2.2. The cross section. Here we define the cross section, Ai, of the "full" billiard flow phase 
space with respect to which we are eventually going to define the billiard maps. The cross section 
itself will be independent of the location of the white disk. For that reason the white disk will not 
enter the discussion at this stage. 

We begin by labeling the non-random walls Ti, . . . , r§ of the domain as shown in Figure 5. The 
boundary of the gray disk forms the four solid walls Tj with an odd index i, and the other four 
with an even index are the transparent walls. The walls are "closed" in that they contain their 
endpoints. As customary, let the position r on a wall be parametrized by arclength and denote by 
ip the clockwise angle relative to the normal vector of the wall pointing into the domain. Facing in 
the direction of the normal vector, the value of r increases from left to right. Then any unit vector 
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Figure 4. The free zone condition. The straight lines are tangent to the gray disk. 
The white disk is not allowed to enter the region with a light gray shading. The 
distance from the center of the square to the latter region is L. 

based on any of the walls and pointing into the domain corresponds, in a one-to-one manner, to 
an element x = (r, (p) of the disjoint union 

8 

A* = Hr,x [-!,§]. 

The mentioned cross section, Ai, of the billiard flow phase space will be obtained by deleting 
a part of Ai as follows. Pick an arbitrary transparent wall Tj, i even, and consider an arbitrary 
point x = (r,ip) G Tj x [— |, |] representing a velocity vector on IV If x corresponds to a crossing 
which does not yield a clean pass (see Figure 2) as described in the Introduction, we delete x from 
the coordinate rectangle Ti x [— |, |]. Such values of x are illustrated in Figure 5. The remaining 
parts of the coordinate rectangle form a region reminiscent of a rhombus, whose closure we denote 
by Aii. We have thus defined 

Mi C Ti x [-|, f], i even. 
On solid walls we include all velocity vectors, thereby setting 

Mi = TiX [-§,!], i odd. 
Finally, the cross section itself is defined as 

8 

M = \\Mi c Ai. 

i=i 

In plain words, each collision with the gray disk marks a return to Ai, in that the post-collision 
velocity vector is an element of Ai. On the other hand, a crossing of a transparent wall is a return 
to M. if and only if it yields a clean pass (see Figure 2). Thus, the three crossings of T§ shown in 
Figure 5 are not returns to Ai. 

It is clear (see Figure 5) that there exists a constant d > such that 

cosy? > d if x = (r, cp) G Ai is on a transparent wall. (1) 

Indeed, the interpretation of the above discussion is that the particle is stopped (x G Ai) at a 
crossing of a transparent wall only when the crossing angle is sufficiently far from being parallel 
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(depending on position) to the transparent wall. In the opposite case the particle is allowed to 
pass through without registering a return. 




Figure 5. The cross section Ai of the billiard flow phase space does not contain 
the shown vectors on r 6 , because they correspond to "immediate" collisions (dots 
on the solid walls) with the gray disk. Deleting all such points from the coordinate 
rectangle leaves us with the white rhombus-like region AA% C Ai (sketched above), 
corresponding to all clean passes through the transparent wall (Figure 2). 



2.3. The billiard maps. Let c denote the position vector of the center of the white disk relative 
to the center of the square. Given any c and 6 > 0, we will throughout the text denote by 
B$(c) the ^-neighborhood of c. We say that c is admissible, if c G B e (0). Having defined the 
cross section Ai C Ai, the billiard map F c is defined simply as the first return map from the 
cross section Ai to itself of the billiard flow $ c corresponding to the white disk centered at c. If 
T c : Ai — > M. stands for the return time function, then 

F c = : M ^ Ai. 

The map F c is called admissible if c is admissible. Notice that, due to the geometry, T c (x) is 
defined and uniformly bounded above for all x G Ai. Also F c (x) is well defined, save for relatively 
few values of x for which the billiard trajectory starting from x meets a corner point, resulting in 
conflicting candidates for the value of F c (x), one on each of the three walls meeting at the corner. 
The conflict is resolved by identifying the three values. Such exceptional values of x are examples 
of singularities, and we will come back to their role shortly. For now, it is harmless to ignore them. 

It will be technically beneficial to view the map F c in a different way. Namely, notice that the 
billiard trajectory between any x G Ai and its image F c (x) G Ai either misses the white disk 
completely or hits it precisely once. Now, let F* be the billiard map which, in addition to Ai, 
counts collisions with the white disk as returns. More precisely, let Ai* be the enlarged cross 
section Ai* = Ai II T* x [— 7r/2, 7r/2], where T* is a parametrization of the boundary of the white 
disk. The parametrization can be fixed in such a way that Ai* is independent of c, and this is 
what we do. Then F* : Ai* — > Ai* is the first return map to Ai* of the billiard flow, defined in a 
fashion similar to F c . Next, let n c : Ai* — > {1, 2} be the smallest number of returns to Ai* which 
yields a return to Ai, i.e., 

n c (x) = minjri > 1 : (F*) n (x) G Ai}, x G Ai* . 
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It is clear that n c (x) = 2 if x G Ai and the billiard trajectory connecting x to F c (x) experiences a 
collision with the white disk. Otherwise n c (x) = 1. With these definitions, 

F c = (F;y°\ M . 

Figure 6 illustrates the situation. 




Figure 6. The action of the billiard map F c corresponding to the white disk being 
centered at c. We have x' = F c (x) and y' = F c (y). On the other hand, z' = F*(z) 
and z" = F*{z') = F c {z). 

We will use the shorthand notation {n c = n} = {x G Ai* : n c (x) = n} with n = 1,2. Notice 
that 

M* = K = 1} U {n c = 2} 

and 

K = 2} = (F*)~ 1 (M* \ M) C M (2) 
hold true by construction, and that 

f;(k = 2})cK = i}. 

2.4. Invariant measure. We may interpret any xGjVfas the velocity vector of a billiard particle. 
Reversing the direction of the velocity to its opposite amounts to mapping x to another point 
X(x) e M. It is easy to check that the map X is an involution on Ai, i.e., it is bijective (taking 
the earlier identifications at the corners into account) and X 2 is the identity map of Ai. 

Picking an arbitrary element y G A4, it is possible to trace the billiard trajectory leading to y 
backwards, by following the trajectory starting from T(y) forwards. By construction, the latter will 
eventually "make a return" to Ai, precisely at F c (X(y)). Reversing the direction of the velocity 
once more, we obtain the element x = X(F c (X(y))) G M. It is clear that F c (x) = y and that x 
is the unique preimage of y. This is to say that the billiard map F c : M. — > M. is bijective, the 
inverse billiard map being given by 

F- l =loF c ol. (3) 

It is also a standard fact that F c = preserves the probability measure 

d/i(r, (p) = M -1 cos ip dr dip 

on Ai, where M = j M cosip dr dp is the normalizing factor, for all admissible c. (This is a 
consequence of the general fact that Hamiltonian flows preserve the Liouville measure.) We point 
out that both of these properties fail generically if the cross section is constructed on the boundaries 
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of moving scatterers; see [26] for examples. (The failure of invariance of p has been overlooked 
in [28].) This is the first of the two fundamental reasons we defined the cross section A4 in 
the specific way that we did, and resorted to using non-moving transparent walls. The second 
fundamental reason is related to guaranteeing uniform hyperbolicity of the maps F c , as explained 
in Section 4.2; see Remark 7 in particular. 

3. Main results 

Given a sequence (c„) n > of admissible centerings, we write 

T n (x) = T n ((c n ) n > ,x) = F Cn _ 1 o • ■ ■ o F Co (x), 

for short, with the convention that J-q is the identity map id^j of Ai 
invariance principle has to do with approximation of the sums 

n— 1 n— 1 

A *(( c i)i>o> x ) = X^ f ( Ci '^(( c ^^°' x ))' n - °> 

i=0 i=0 

suitably scaled, by Brownian paths in M. d , given an Revalued observable f together with a ran- 
dom sequence of admissible centerings c„ of the white disk. The strength of the statement is in 
that the approximation converges almost surely (contrasting the weak invariance principle which 
only guarantees that the finite-dimensional distributions converge.) Our result will concern, pos- 
sibly non- stationary, random sequences (c n ) n > , where the random variables c n can be weakly 
dependent. Before stating the theorem, we introduce the assumptions on the sequences (c n ) n > . 

First, we recall a notion of weak dependence from probability theory. Suppose a probability space 
and two sub-sigma-algebras 21 and 33 are given. The so-called maximal correlation coefficient of 
21 and 33 is given by 

p(Sl,») = sup{|Corr (/,(/) | : / G L 2 (2l;IR), g G L 2 (Q3;IR)}. 

It can be shown [4, 29] that 

p(2l,^) = sup( |E(/ ^~^ ( { )E( ^ )l : /GL 2 (2l;C),seL 2 (<B;C)l. (5) 

(Here L 2 (2t;lR) is the space of real-valued, 2t-measurable, square-integrable functions, etc.) Note 
that the supremum in (5) is taken over complex-valued functions and that < p(2l, 53) < 1. 
Consider next a random sequence (X n ) n > and the probability space of its trajectories. For n > 
m > 0, let #m be the sub-sigma-algebra generated by the variables X m , . . . ,X n . The so-called 
rho-mixing coefficients are 

p(*)=supp(SS,3ft*), k>o. 

n>0 

The sequence (X n ) n > , or its probability distribution, is called rho-mixing, if lim^oc p{k) = 0. 
Obviously, a sequence of independent random variables is rho-mixing with p(k) = for k > 1. 
More generally, the same is true of M-dependent sequences, for any M G N. A stationary Markov 
chain is rho-mixing if it has the L 2 -spectral-gap property, and in this case p(k) tends to zero 
exponentially [22]. Moreover, if a Markov chain, stationary or not, satisfies p(ko) < 1 for some 
k > 1, then p{k) tends to zero exponentially [4,22]. 

Given e > and a probability distribution P of the random sequence (c n ) n > on (B e (0)) n 
endowed with the product Borel sigma-algebra our first assumption is the following. 



x e M, (4) 
. The K -valued almost sure 
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(Al) P is exponentially rho-mixing. That is, there exist c > and B > such that 

p(k) < Be~ ck , k>0. 

Assumption (Al) is sufficient in the stationary case. In order to deal with non- stationary 
sequences, we impose a further condition which holds automatically in the stationary context. It 
is related to the concept of asymptotic mean stationarity (see, e.g., [9] and the proof of Lemma 2 
below), which means that the measures F k = | ^^(cr- 7 )* P, k > 1, obtained by averaging over 
shifted sequences, tend to a measure P which is (necessarily) stationary. Here and below, a is the 
left shift on (B £ (0)f. 

Given a measurable function g : (B £ (0)) N — > R, let us introduce the shorthand notation 

fe-i 

(g) k = -Y,E(gooS), k>l, (6) 

j=0 

where E is the expectation relative to P. Our second assumption is the following. 

(A2) For any bounded measurable function g, there exists a number (g)oo G K such that lim^oo (g) k = 
(g)oo- There exist sequences (C m ) m >o and (rk)k>o of positive numbers with r\ — >■ such that, for 
any m > and any bounded ^-measurable function g m , 

\{g m ° <y l )k - (sw)oo| < C m r k Halloo, k>l,£>0. (7) 
Finally, for any (3 > 0, there exists a constant Dp > such that 

/3~ ll0gn J g n 

/ 1 C m r n _ m e ^ m < (8) 
^— ' n 

holds whenever n > logra (so that r n _/3-ii ogn on the left side makes sense). 

We make the remark about assumption (A2) that if {g m o a e )k and (g m )k converge as k — > oo, 
they obviously converge to the same limit (Lemma 20) for all i > 0, but (7) is needed for a uniform 
rate of convergence and (8) for the rate to be sufficiently fast. 

We are in position to state our main result. As usual, also our billiard maps have singularities 
consisting of a finitely many smooth curves in Ai. For a (non-random) c, denote by S c the 
singularity set of F c ; see Section 4.1. For later convenience, we switch to writing uj = (u n ) n > for a 
realization of the random sequence (c n )„> . That is, u n denotes the realized position of the center 
of the white disk at time n. 

Theorem 1. There exists a number e > such that the following assumptions have the following 
consequences. 

- Assumptions: Let P be a probability distribution on (B e (0)) N satisfying (Al) and (A2). Write 
P = P ® a and E( ■ ) = /( ■ ) dP. Let Q C B e (0) be a measurable set with P((fi ) N ) = 1. Let 
f : flo x M. — > lR d (any d > 1) be a bounded measurable function such that, for all c e Qq, f(c,x) 
is defined for \i-almost- every x G M. and f f (c, x) d/i(x) = 0. Let the maps f(c, •), c G fi 0; be 
uniformly piecewise Holder continuous in the sense that there exist 7 > and Cf > for which 

\f(c,x)-f(c,y)\<C t d(x, y y 

for allx,y in the same component of M. \ S C , for all c G Qq. Denote Ak(uj,x) = f (uk, J^k^, #)) . 
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Consequences: 

(1) The formula 

A— J 

1=0 



S '= lim f (A ® A )(^o;,a;)dP(a;, 

fe ^OO f£ I 

1=0 J 

X ^ fe— 1 „ 

+ V lim - V / (A ® A m + A m <g> A ) (> £ w, ar) dP(w, ar) 

fe — ^oo fZ I 



m=l 1=0 

yields a well-defined, symmetric, semi-positive-definite, d x d matrix E 2 . 

(2) The matrix E 2 is the limiting covariance of X]fe=o Afe- That is, 

- / n— 1 n— 1 \ 

lim -E VA fc ® VA fc = E 2 . 

n— >oo 77, \ ' ^— ' / 

\fc=0 fc=0 / 

(3) The random variables X]fc=o Afe converge in distribution, as n — > oo, to a centered IR d - 
valued normal random variable with covariance S 2 . 

(4) Given any A > \, there exists a probability space together with two W 1 -valued processes 
(A* )„> and (B„) n > on it, for which the following statements are true: 

(a) (A n )„> and (A*) n > have the same distribution. 

(b) The random variables B„, n > 0, are independent, centered, and normally distributed 
with covariance S 2 . 

(c) Almost surely, | ^^=0 A t ~ Efc=J B fcl = 2 

Item (3) of the theorem is called the averaged central limit theorem and item (4) the vector-valued 
almost sure invariance principle with limiting covariance S 2 and error exponent A. The "almost 
surely" in item (c) refers to the probability space on which the processes (A*) n > and (B n ) n > are 
defined. Since Y^k=o nas ^ ne interpretation of the location of an IR d - valued Brownian motion 
at time n, the result is indeed about Brownian approximations of the sequences (^^=0 Afe)«>o- 
Accordingly, [3,20,27] discuss several implications of the almost sure invariance principle, and [10] 
shows that the almost sure central limit theorem is among them. 

Note that the centerings c n in the theorem could for example be drawn independently, from the 
same distribution; the distribution could be uniform in the disk of radius e, which in a sense corre- 
sponds to maximal randomness, or it could be a delta-distribution corresponding to a completely 
fixed white disk as in Sinai billiards. 

Regarding the regularity condition imposed on the observable, the inter-return flight time and 
vector displacement of the particle lifted to the plane both satisfy it [6]. 

The expression for the limiting covariance in (9) is quite remarkable in that it applies to the 
non- stationary case. In addition to the expectations with respect to the distribution P = P <g> u, 
the terms in the infinite series entail averaging over the trajectories (a e u)e>o of the sequences u 
under the left shift a. In the stationary case the averaging is redundant and the expression of the 
covariance simplifies to the well-known series E(A <8> A ) + X^m=i E(Ao ® A m + A m <g> A ). The 
matrix S 2 can be degenerate for a given observable f . To characterize such a situation, we prove 
the next lemma in Section 5.6, following the classical works of Robinson [21] and Leonov [13]. 



As is customary, given a sequence (a n )„>o, the notation a n = o(n x ) means that lim„ 
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Lemma 2. Let f and E 2 be as in Theorem 1. Denote by $ the skew-product map defined as 
x) = (era;, F uo (x)). The matrix S 2 is degenerate if and only if there exists a (constant) 
vector v G M d and a measurable function g : (fio) N x M. — > K such that, for all j > 0, the map 
(u,x) I—)- g(a^uj,x) belongs to L 2 (P), lim^oo |X^j=oJ IS'O 7 *'^) x )\ 2 dP(w, x) < oo, and t/ie identity 
v T f = g — goQ holds except on a measurable set E C (fio) N x A4 which is asymptotically negligible 
in the sense that lim^oo | X}j=o P( a = 0- 

Remark 3. Theorem 1 does not depend on the assumption that the scatterers are disk-shaped. The 
same proof continues to work for more general, strictly convex, shapes of the fixed "gray" scatterer 
and the moving "white" scatterer having at least C 3 -smooth boundaries, because [26] applies to 
such geometric objects just the same. Of course, it is still necessary that the free path length 
between the solid scatterers is uniformly bounded from above and below by positive numbers. (This 
was guaranteed by the no-overlap and finite horizon conditions above.) In order not to have to 
upgrade the proof, the billiard trajectory should also not hit the white scatterer more than once ( or, 
a bit more generally, a uniformly bounded number of times) between successive returns to the cross 
section AA. (This was guaranteed by the free zone condition.) The last thing one has to make sure 
of is the preservation of uniform hyperbolicity despite the transparent walls by a suitable choice 
of AA; between successive returns involving one to a transparent wall, the flight time and the angle 
of the velocity at the transparent wall should be uniformly bounded. (This followed from the clean 
pass condition; see {!)■) Finally, besides just translations, rotations as well as changes to the shape 
of the white scatterer can be allowed, provided they are uniformly sufficiently small, because [26] 
applies equally to that scenario. We leave the formulation to the reader. 

The next result concerning fixed scatterer configurations (Sinai billiards) and random observables 
can be deduced from the very same proof as Theorem 1, simply switching to the familiar cross 
section Af corresponding to returns of the billiard flow to the solid boundaries of the scatterers. 
Since the cross section to be used is clear from the beginning, we formulate the theorem in a rather 
general geometry (not assuming circular boundaries). We also assume a very modest amount 
of regularity of the observables, namely dynamical Holder continuity in the stable and unstable 
directions separately. (In particular, "ordinary" Holder continuity suffices. See Section 5.1 and 
Remark 10 there for a discussion of dynamical Holder continuity.) 

Theorem 4. Consider a Sinai billiard on the two-dimensional torus with strictly convex scatterers 
having C 3 -smooth boundaries. Assume the free path length is bounded from above and below by 
positive numbers. Let F : Af — >■ Af denote the usual billiard map; see the paragraph above. Let Q 
be an index set, d > a fixed integer, and f : Qq x Af — > IR d a bounded measurable map with 
the property that (each vector component of) the map f(wo, • ) : Af — > M. d is dynamically Holder 
continuous on homogeneous local stable and unstable manifolds, with uniform parameters for all 
ujq G VLq. Let ¥ be a probability distribution on the space (fio) N of sequences u satisfying (Al) 
and (A2). Denote A k (u, x) = f{u k , F k (x)), k > 0. Also write P = P <g> fi and E( • ) = /(-) dP - 
Then the consequences listed in Theorem 1 and Lemma 2 hold. 

The author is unaware of other limit theorems for Sinai billiards with time-dependent (here 
random) observables. Turning to the special case of a fixed observable (that is, Qq contains a 
single element), the above theorem enlarges the class of observables for which a weak or almost 
sure invariance principle is known [5,15,16]; the central limit theorem for the same class was proved 
in [24]. A virtue of the result is that, besides being quite large, the class of observables in question 
is convenient to work with when it comes to the study of statistical properties; see discussion 
below. In addition, Theorem 4 greatly improves the best known error exponent for scalar-valued 
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observables (from | + e to \ + e) as well as for vector-valued observables in all dimensions d > 2 
(from ||±| + e again to \ + e) [16] . 

Giving a separate proof of Theorem 4 would be pointless, because that would amount to repeat- 
ing the exact proof of Theorem 1 with the simplification that the map and the function spaces 
would no longer depend on time. Accordingly, it is also possible to formulate a more elaborate 
(albeit general) version of Theorem 1, in which the observables belong to time- dependent spaces 
of dynamically Holder continuous functions (see Section 5.1); after all, the proof itself is written 
precisely for such spaces. But for the sake of a more intelligible statement, and to stress what is 
really relevant and new, we offer Theorem 1 in its current form. 

The proof of Theorem 1 benefits from a condition formulated recently by Gouezel under which 
the vector-valued almost sure invariance principle holds for a sequence of random variables. We 
will recall the condition in Section 5.4. Gouezel's proof relies on a sophisticated variation of a 
classical blocking argument, in which the marginals of the random process are grouped into blocks 
separated by a time gap of length k. His condition quantifies a sufficient rate for two blocks to 
become independent, as their separation k increases, for the vector- valued almost sure invariance 
principle to hold. 

The strategy of verifying Gouezel's condition in the proof of Theorem 1 involves a mixture of 
ideas related to [24,26]. Viewing the model under study in a way amenable to analysis, including 
how we defined the cross section Ai and the billiard maps F c on it, is also very much part of the 
method. (See Section 4 for more on the importance of setting up the problem carefully.) In the first 
reference, [24], the author enlarged the function spaces considered in [6] by relaxing the regularity 
assumptions. To include less regular observables in the analysis of the statistical properties of 
billiards turned out to be fruitful for proving limit theorems. It was shown that, for observables 
in these enlarged spaces, a pair correlation bound alone implies the central limit theorem. In the 
present paper we take advantage of similar large spaces, but this time in a time- dependent setup, 
and prove bounds on correlation functions which imply Gouezel's condition for the invariance 
principle. The theory of time-dependent billiards has been developed in [26] to the stage that the 
existence and uniformity of hyperbolic structures, for example, can be taken for granted in the 
analysis to follow, given the deliberate way in which we defined the cross section Ai and the billiard 
maps F c . Whereas the general paper [26] did not focus on random compositions of billiard maps 
let alone on their limit theorems, the latter form exactly the goal we are now shooting for. 

4. Preliminaries 

In this section we introduce some concepts and facts that enter the proof of Theorem 1. We are 
not going to present proofs or new results, save for Lemma 6. Due to the slightly unusual choice 
of Ai , we need to verify the uniform hyperbolicity of the maps F c — in the form of Lemma 6 - 
by hand. After this has been achieved, we are immediately in position to apply the machinery of 
time-dependent billiards developed in [26]. 

4.1. Singularities. Like in the case of classical billiards, the maps F c have singularities. The 
qualitative behavior of a billiard map changes in arbitrarily small neighborhoods of its singularities. 
Given a centering c, we call 

S c = F-\dM) U (F*)-\dM* \ dM) 

the singularity set of F c . The first member of the union corresponds to those values of x G Ai for 
which the billiard trajectory from x to F c (x) either hits an endpoint of a wall (corner); meets the 
gray disk or a transparent wall tangentially; or crosses a transparent wall in such a way that it is a 



14 



MIKKO STENLUND 



return to Ai, but there are arbitrarily small perturbations of the trajectory which are not returns 
to Ai. The second member of the union corresponds to those x £ Ai for which the trajectory 
meets the white disk tangentially. We point out that (2) yields (F*) _1 (9A^* \ dAi) C Ai, which 
is to say that S c as defined above is a subset of Ai as it should. 
The reader will notice that 

S c = {F* c \ M y\dM*) U (F^m)- 1 (f:\m.\m)- 1 (dM), 

which is a convenient way of viewing the singularity set. The interpretation is that the second 
member of the union corresponds to "secondary" singularities — singularities as first seen from 
the white disk and then pulled back to Ai. 

The singularity set S c C Ai consists of piecewise smooth curves with uniformly bounded negative 
slopes. (In fact, they are stable curves, i.e., curves S whose tangent spaces T X S, x £ S, are 
contained in the respective stable cones C x as defined below.) The smooth curves, or branches, 
either terminate on other branches in the interior of Ai , or they extend all the way to the boundary 
of Ai, and every branch is a part of a path of branches reaching from one part of the boundary 
to another, as is well known. When c varies, the singularity set S c varies in a continuous manner: 
Its topology can change in that new branches can appear and existing ones can disappear, and 
the branches generally undergo deformations, but S c is contained in an arbitrarily small tubular 
neighborhood of 5 C provided c is in a sufficiently small neighborhood of c. Moreover, the number 
of branches is uniformly bounded from above. See [26] for more on time-dependent singularity 
sets. In particular, it is possible to prove [26] that 

Lemma 5. Given an admissible c and a compact subset E C Ai \ Sc, there exists 5 > such that 
the maps (c, sc) h> F c (x) and (c,x) t— > D X F C are uniformly continuous on E x Bg(c). 

The map (c, sc) i-> F c (x) is thus continuous on the set {(c, sc) £ B £ (0) x Ai : sc £ Ai \ S c }. 

We can extend the concept of singularities for sequences of billiard maps in a natural way. 
Suppose, then, that c and c x are given centerings. The set 

<Sc ,ci = «S Co U F Co («S Cl ) 

is the singularity set of the composition F C1 o F CQ . More generally, given a (finite or infinite) 
sequence (cj)^, 1 , we can define inductively, for 1 < n < N , the singularity set of the composition 
F Cn _ ± o • • • o F co as 

<S C0 ,...,c n ^i — <S CQ U F CQ («S clr .. jCn _ 1 ). 
The importance of these sets lies in the fact that the actual billiard trajectory from sc to F Cn _ 1 o 

• • • o F Co (x) meets a singularity at some point if and only if x £ <Sc ,...,c n _i- 

Going backward in time, similar singularity sets can be defined for the compositions (-Fc„_i ° 

• • • ° -Pen) -1 - We denote them S~ „ , . 

4.2. Uniform hyperbolicity. Let sc = (r, <p) £ Ai be arbitrary. We denote by /t gra y and ft w hitc 
the curvature of the boundary of the gray and the white disk, respectively. If sc is on a solid wall 
(of the gray disk), define 

cos 

a(x) = K gray and b(x) = K gray H . (10) 

Tmin 

On the other hand, if sc is on a transparent wall, define 

a(sc) = — and 6(sc) = — , (11) 
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where K min = min(K gray , K whitc ). Note that b(x) is strictly positive by (1). In fact, there exists 
constants a min and 6 max such that 

< a min <a(x) < b(x) < 6 max < oo 

for all x G M. 

Lemma 6. For any admissible c, the unstable cones 

CI = {(dr, dip) G T X M : a(x) < d^/dr < b(x)}, x G M, 

are invariant for F c and the stable cones 

C s x = {(dr,d<^) G T X M : -a(x) < d^/dr < -b{x)}, x G M, 

are invariant for F~ l . There exists constants C > and A > 1 such that, given n > 1 admissible 
centerings c , . . . , c n „ 1; the bound 

\\D x (F Cn _ 1 o...oF CQ )v\\>CA n \\v\\ 

holds for all v G C x and for all x G M. \ «S COj ... jCn _ 1 , and the bound 

P^Fe^o-.-oFeJ-SlI >CA n ||t;|| 

holds for all v G Ci and for all x G M. \ S' 1 _ 

The lemma states that the maps are uniformly hyperbolic with common families of unstable 
and stable cones and with a uniform minimal expansion rate in the cones. We are therefore in a 
setting to which the time-dependent analysis that was carried out in [26] applies directly. 

Proof of Lemma 6. Consider a point x\ = (ri, ipx) G M* and set xq = (r , <^o) = (F* G M* 
which corresponds to a straight flight from x to X\ (possibly passing through a transparent wall 
in a way which does not constitute a return to M). We have [6] 

V. = Kl + -^. (12) 

Here tq is the path length from xq to xi, /tj = tz(xi) is the curvature of the wall at Xj, Vo is the 
slope of a tangent vector (dr, dip) at x and Vi is the slope of the corresponding tangent vector at 
x% obtained via the tangent map. Observe that, by our conventions of defining the cross section, 
we have the uniform bounds 

< r min < r < r max < oo. (13) 

First, we see that the "positive" unstable cones {(dr, dip) G T X A4* : < dip/dr < oo}, x G M*, 
are invariant for F*. In fact, the image of a "positive" tangent vector at xq satisfies 

Kl < Vi < Ki H . 

We wish to improve the lower bound on Vi in the case k± = 0, i.e., when x\ is on a transparent 
wall and in M.. In this case x , defined above, is necessarily on a solid wall. Thus, k > n min > 
and 

cos^o < 1 



V + 

Recalling (1), we have cos^i > d, so that, by (12) and (13), 

coswi d 
Vi = — — > 

y J- , COS (Z)n 



cosyp - ^_ 
'0 -T Vo+KQ /max T 
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Let us extend the definition of from x G Ai to all x G Ai* by defining a(x) and as in (10) 
for x G Ai* \ Ai. By the above analysis, we conclude that F* preserves the obtained cones C", 



Finally, suppose x is a general element of At and set x\ = F c (x) G Ai. Then the billiard 
trajectory from x to Xi can be broken into straight legs, the last one of which is from Xo = 
G Ai* to x\. Since each of the legs preserves the cones C", x G A4*, the map F c 
preserves the cones C", x G A4. 

Coming to the stable cones, we resort to a time-reversal argument. Recall the expression (3) of 
the inverse billiard map F~ x in terms of the involution X. As can be checked, the derivative D X X 
equals either (o -i) or ( V ?) depending on whether i is on a solid or a transparent wall, re- 
spectively. Using the obvious identity D X I(C%.) = Cx{x) anc ^ the above fact that F c preserves the 
unstable cones, one confirms that F~ l indeed preserves the stable cones. 

Next, we prove the uniform expansion property of the compositions F Cn _ 1 o • • • o F Co on unstable 
vectors. The complementary statement concerning stable vectors and compositions of inverse 
maps is similar. To that end, fix a y G Ai and denote = F Ck _ 1 o ■■■ o F CQ (y ) G Ai for 
1 < k < n. We assume y ^ <5 COi ... )Cn _ 1 . Also fix an unstable vector dy G Cy Q and denote 
(h)k = -Dj /0 (F Cfe _ 1 o • • • o F Co ) dyo for 1 < k < n. The billiard trajectory from y to y n can be broken 
into at most 2n straight legs: Given 1 < k < n, yu = F Ck l {yk-i) equals either F* (F* (yk-i)) or 
F* k (yk-i) depending, respectively, on whether the trajectory from y k _ x to y^ involves a collision 
with the white disk or not. Thus, we obtain collision points Xj G Ai*, < j < m, where 
n <m < 2n. Of course, {xj)^ =0 has (yk)k=o as an ordered subsequence, with x = y Q and x m = y n . 
Likewise, we get a sequence of tangent vectors dxj, < j < m, so that dxo = dyo and dx m = dy n . 

The billiard maps F* k of the extended cross section At* satisfy, by well-known [6] formulas, 



Here || • || p is the so-called p-metric of the tangent space, Vj refers to the slope of dxj, Kj to the 
curvature of the boundary at Xj = (r^, (fj), and Tj to the length of the straight leg from Xj to Xj + ±. 



x G Ai*. 



ll da; i+i||p 
lld^jllp 



1 + Tj 



Vj + Kj 



COS (fj 




dx\\ = ^Vi^v\ 




||dy n || ||dx m || ||dx m || p ||dxi|| p cosv?o a/TTV^ 




(14) 



Therefore 
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with 



which proves the lemma. 



C = A 



-i 



T, 



min^mm 



1 + b 2 



1 + a 2 

1 T "mm' 



□ 



Remark 7. Notice in the proof that if the image point x% was on a transparent wall, K\ = 0, and 
if in that case it was possible for cos 921 to be arbitrarily small in contrast to (1), then (12) would 
result in arbitrarily small values of the slope V\ of the image tangent vector. Subsequently, we 
would not obtain a positive lower bound a m ; n in (11). Even worse, if it was possible for cos<^o to 
be arbitrarily small when Xq is on a transparent wall, we would not obtain a positive lower bound 
in (14), nor the uniform expansion rate A > 1 in Lemma 6. These facts constitute the second 
fundamental reason, alluded to in Section 2, why we defined the cross section Ai in the specific 
way, ruling out small values of the cosine on transparent walls. 

4.3. Homogeneous local stable manifolds. In order to control distortion effects of dispersing 
billiard maps, it is traditional to introduce so-called homogeneity strips. In our case, 



U k = {(r,<p) G M* : tt/2 - AT 2 < up < tt/2 - (k + 1) 
= {(r, if) G Ai* : -tt/2 + (k + 1)~ 2 < if < -tt/2 



2 } 

hk~ 2 } 



for all integers k > ho, where ko is a sufficiently large uniform constant. We also set 

H = {(r, if) G Ai* : -tt/2 + A; 2 < if < tt/2 - kf} . 

Without going into any detail, if F*(x) and F*(y) belong to the same homogeneity strip, then 
roughly speaking the derivatives D X F* and D y F* are comparable. By (1), we may assume that k 
is so large that each component Aii with i even (transparent wall) involves just one such strip, 
namely Ho- 

Suppose that (c n ) n > is a sequence of admissible centerings. We say that two points x and y 
in M* are separated if they lie either in different connected components of Ai* or in different 
homogeneity strips of the same component. Recall the notation in (4). Given a pair (x,y) G 
M. x A4, the future future separation time s + (x,y) = s + ((c n ) n > ; x, y) of x and y is defined as 
follows. If x and y are separated, s + (x,y) = 0. Otherwise s + (x,y) is the smallest integer n > 1 
such that either F* n ^{F n -i{x)) and F* ji _ i (J r n _i(?/)), or J- n (x) and F n (y), are separated. (The first 
alternative takes into account the possibility that the trajectories land on different homogeneity 
strips on the white disk, or that only one of them hits the white disk, before making a return to 
M.) 

Let (c ra ) n > be an infinite sequence of admissible centerings. A connected smooth curve W is 
called a local stable manifold, if the following hold for every n > 0: (i) F n {W) is connected and 
homogeneous, (ii) T x (F n (W)) C C s x for every x G F n (W), and (iii) similar statements hold for 
F Cn ((F n (W))). (Item (iii) takes into account the possibility that the curve hits the white scatterer 
before returning to Ai.) The local stable manifolds are absolutely continuous with uniform bounds 
on the Jacobian, the key factor behind it being the uniform hyperbolicity discussed before. By the 
same token, they contract at a uniform exponential rate (Lemma 6). Notice that, by definition, 
two points on the same homogeneous local stable manifold never separate, so that their future 
separation time is infinite. It turns out that this is a defining property of homogeneous local stable 
manifolds. See [26] for more on local stable manifolds in the time-dependent setting. 
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4.4. Exponential loss of memory. Given any admissible centering c, classical theory of billiards 

shows that the map F c , corresponding to the fixed center c, is ergodic and mixing with respect to 
the invariant measure p. 

For the compositions (4) of such maps with time- dependent centerings c n , the analysis carried 
out in [26] is directly applicable. The main result of that paper can be stated in the present context 
as follows: 

Theorem 8. There exists e > such that the following holds. Let p 1 be a probability measure on 
A4, with a strictly positive, ^-Holder continuous density p 1 with respect to the measure p. Given 
7 > 0, there exist < # 7 < 1 and C 7 > such that 



for all sequences (c n )^l of admissible centerings and all ^-Holder continuous f : Ai — > R. The 
constant C 7 = C 7 (p 1 ) depends on the density p 1 through the Holder constant of log p l , while # 7 
does not depend on p 1 . 

Noting that f M f o T n dp 1 = f M f d(J r n )*/i 1 , the theorem states that the push-forward measure 
{J r n)*p l converges at a uniform exponential rate to p, in a properly understood weak sense. In 
particular, recalling that mixing can be formulated as correlation decay (see Theorem 11), when 
applied to constant sequences c n = c, n > 0, the above theorem yields a uniform exponential 
mixing rate for Holder continuous observables, for all fixed admissible maps F c . 

Theorem 8 has a pre-eminent role in the proof of Theorem 1. In truth, we will need a slightly 
refined version of it (Lemma 12) which allows for the measure p 1 to be supported on a single 
"unstable" curve. Nevertheless, it seems appropriate to introduce the above formulation as a key 
ingredient of the theory. 



The proof of Theorem 1 — as is usually the case with billiard proofs — is rather technical. 
In order to keep the paper in manageable proportions, we will need to assume that the reader 
has a working knowledge of billiard techniques. Nevertheless, the flow of the argument should be 
accessible to a nonspecialist. For background on billiards, we refer to the excellent textbook [6] 
and, on the time-dependent theory, to the paper [26]. 

The proof will proceed as follows. In Section 5.1 we introduce some notations and define the 
spaces of observables. In Section 5.2 we obtain uniform, exponentially decaying, bounds on pair 
correlation functions for observables in those spaces. In Section 5.3 we extend the result to multiple 
correlation functions. Using the multiple correlation bounds, we can control the characteristic 
function of the process, and thus verify Gouezel's condition for the vector-valued almost sure 
invariance principle. This is done in Section 5.4. In the non-stationary case, we also need a bound 
on the converge rate of the covariance, which is obtained in Section 5.5. Finally, Lemma 2 on the 
(non-) degeneracy of the limiting covariance is proved in Section 5.6. 

5.1. Some definitions. Here we will denote by 



the space of bi-infinite sequences of admissible centerings. Elements of f2 are denoted by oo = 
(wi)ieZj where each Ui is thus an admissible centering inside the disk B £ {Q). The usual left shift is 
denoted by a, i.e., 




n > 



5. The proof 



Q = (B E (0)) 



a : Q — > Q : (au))i = u i+ i Vi G Z. 
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For technical reasons, we will work with bi-infinite sequences, although the sequences (F Ui )i> Q is 
what we are eventually interested in. 

Given a sequence uj G Q, we will use notation such as 

F n (w,-) = F Un _ 1 O'--oF U0 (-) and T n+m>n+1 (u, ■ ) = T m (a n u, ■ ), 

as well as 

K V> • ) = (^n-i ° • • • ° F^T\ ■ ) and 7^ TO)n+1 (a;, • ) = T^{a n u, ■ ) 

for n,m > 0. We use the convention that J-^n+i = idx throughout. The notation is consistent 
with that, given u, the maps F Wn and .F„ +m) „ + i(u;, • ) describe the transformations from time n to 
time n + 1 and to time n + m, respectively. 
Let 

be the measurable partition of Ai consisting of the homogeneous local unstable manifolds corre- 
sponding to the past (. . . , U-2, uj-i) and 

V-{u) = (W-(u) a ) aeA _ (u) 

the measurable partition of M. consisting of the homogeneous local stable manifolds corresponding 
to the future (oj , ui, . . .). Also set 

V { ±\u) =V±(a n u), neZ. 

We often write p£ n) (w) = (w£l(u)) aeA ^ (w) , where 

At\u) = A ± {a n u) and W^u) = W ± , a (a n u), 

or simply 

when there is no danger of confusion about the sequence u. Note that, for any n G 

is a refinement of the partition F u]n (V^\oo)) , obtained by cutting the elements of the latter into 

maximal homogeneous components. Similarly, "P_ (w) is a refinement of F" 1 ^!"^^)). 

Recall that, given a one sided sequence (co> , ^i, • • • ) and two points x,y G Ai, we have defined 
the future separation time of x and ?/ in Section 4.3. The definition obviously extends to two-sided 
sequences u G fl and depends only on Ui, i > 0. We write s + (lu;x, y) for the future separation 
time. Similarly, the past separation time S-(u;x,y) is defined going backward in time, via the 
inverse maps F~} , F~} 2 o F~* , . . . , and it depends only on Ui, i < 0. 

Given a sequence u; and a number $ G (0, 1), we say that a function / : Ai — > C is dynamically 
Holder continuous on the homogeneous local unstable manifolds associated to u with rate if 
there exists a real constant if > such that 

\f(x)-f(y)\<K# s +^ 

holds for all x,y G VK +)Q and all a G ^4 + . The smallest possible such constant K is denoted 
K + (u;$, f). The class of such functions is a vector space, which we denote H + (co;i&). Similarly, 
a function g : Ai — > C is called dynamically Holder continuous on the homogeneous local stable 
manifolds associated to u with rate if there exists a real constant K > such that 

|<7(x)-s(y)| <^-(^) 
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holds for all x,y G W_ )Q and all a G The smallest possible such constant -K" is denoted 

K_(oj] /). The class of such functions is a vector space, which we denote 7i-(uj; 

The above function spaces generalize similar spaces of [24] to the time-dependent setting, whereas 
the spaces of [24] generalize those of [5]. In particular, ordinary Holder continuous functions are 
dynamically Holder continuous: 

Lemma 9. There exists a uniform constant C > such that the following holds. Given Uo G B £ (0), 
7 > and Af > 0, suppose f : Ai — > C is a bounded function that satisfies 

\f(x)-f(y)\<A f d(x, y y 

for all x and y belonging to the same homogeneous component of M. \ Sa, - (Here d( ■ , • ) is the 
Euclidean distance in M..) Then f G H + (u>; A -7 ) fl "H_(w; A -7 ) with 

K + {u;A~\f) <max(A / C' 7 ,2||/|| 00 )A 7 and K^(u; A~ 7 , /) < A f & 'A 7 

for all u £ Q with ojq = uq. 

Proof. Let D be the maximal diameter of the homogeneous components of M.. Suppose u> satisfies 
u = cjq. For any x,y G W +>a and a G A + , CA n d(x,y) < rf(J r n (co', x), J- n (co, y)) < D if < 
n < s + (u; x,y). This follows from the uniform hyperbolicity (Lemma 6). We then arrive at 
d(x,y) < CAA~ s+ ( ul ' x ' y } with the uniform constant C = DC^ 1 . We either have s + (u;x, y) > 1 or 
s + {oj; x, y) = 1. In the first case, x and y belong to the same homogeneous component of M. \ S^ , 
so the assumed Holder continuity of / yields — f(y)\ < A/((7AA~ s+ (' J;a: ' 3/ )) 7 . In the second 

case, \f(x) - f(y)\ < 2^]^ = 2||/|| 00 A 7 (A- 7 ) s +( t ^. The claim concerning H+{uj; A" 7 ) follows. 
Assume now that x,y G W- ;Ce and a G instead. By similar reasoning, d(x,y) < CAA' s -^' x ' y \ 
Forward images of homogeneous local stable manifold do not, by construction, meet singularity 
curves, so that W- >a belongs to a single homogeneous component of M. \Sq q . Thus, \f(x) — f(y)\ < 
AfiCAA-'-^'tiyr. □ 

Remark 10. All of the above applies when the billiard map does not depend on time (fixed scat- 
terers). In that case the partitions V± are defined by the given scatterer configuration instead of a 
sequence, and the same is true of the separation times s± as well as the function classes T-L±. 

5.2. Pair correlation bounds. We say that u G VL has a saturating past, if there exists an i Q < 
such that bJi = Ui for all i < io- Denote by Q sat C Q the set of all sequences in Q with a saturating 
past. Notice that 

°"(^sat) C ^sat) 

that is, the left shift of a sequence with a saturating past has a saturating past. The class fl sat will 
be needed in place of Q because of Lemma 13; see the remark after it. 

The next result gives an exponential upper bound on pair correlation functions for observables 
of the type introduced above. 

Theorem 11. For a sufficiently small e > 0, there exist uniform constants C p > 0, C > 0, 
C > 0, and x > such that, for every $± G (0, 1), every oj G fi sat; every n > 0, and every pair 
(f,g) eH+{u;;0 + ) *H-{a n u>; ■&..), 

j f ■ 9 o 7" n (w, • ) d/i - J fd/i J gd/i 

< WfW^B^u-^g) W-))" /4 + 2C p |mUM|ooe- n/4 * + 2K + (u;$ + J)\\g\\ 00 $ r f-\ 
Here B_(a n uj;$_,g) = CWgW^ + K^(a n uj; 0_, g) $Z l and 0(0_) = max(C, tfl /2 ) G (0,1). 
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Recalling Lemma 9, we point out that the previous applies to the special case of "ordinary" 
Holder continuous observables and yields an exponential rate of pair correlation decay for them. 

As we will see later on, Theorem 1 1 is surprisingly strong in that it implies a bound on multiple 
correlation functions, which is seemingly much more general a result. To understand how this is 
possible, one has to appreciate two aspects of Theorem 11. First, it involves carefully designed and 
sufficiently large classes of observables. Second, the bound is completely explicit in regard to its 
dependence on the observables / and g. In conjunction with Lemma 15 below, Theorem 11 thus 
becomes a very useful tool for proving limit results. 



Proof of Theorem 11. Denoting 



we have 



9 = 9- / 9 dp, 



J f • 9 ° ^n+m(w, • ) d/j, - J fdfi J gdfi = J f o T n ■ ) ■ g o T n+m>n+1 (uj, ■ ) d/x 

for any fixed u. As a; is fixed for the rest of the proof, we will often omit it from the notation. 
Using the partition V+\ the measure \x disintegrates into a probability measure on A { + ] and 
a system of conditional probability measures on the partition elements W+ a such that 

hdn= f [ hdv^d\M(a) 

holds for any Borel measurable function h : M. — > C. We can thus write the above identity as 

J f ■ go T n + m dp - J fdfi J g dfi = h + I 2 

where 



goF n+m>n+1 dvM d\ (n \a) 



and 



h = [ „ tfoF- l du^ ■ (go F nW1 di/W dAW(ot). 

Let us bound I\ first. Because each is a probability measure, 

|ii|<sup sup (f o F~ l {x) -fo F^iy)) ■ / \g\dp. 



For any a G A™ ' there exists a /3 G A+ such that for any pair of points x, y G W+n we have 
J--i(x), F-\y) G W+,p. Therefore, 

Since / G H+(u)\ "&+), we obtain the bound 

|/ o F~\x) - / o T~\y)\ < K + (co; + , fW , 
for all x, y G W + ^ a , for all a G .4.+ . Now sup |#| < 2 sup \g\ yields 

\I 1 \<2sup\g\K{u;# + ,f)#l. (15) 



(n) 
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Coming to J 2 , notice that 



\h\ < sup l/l • / / 00 



n+m,n+l 



dA (n) i 



a). 



(16) 



The following result is at the heart of bounding I 2 . Recall that e is the upper bound on the 
displacement of the center of the white disk from the center of the unit square. 

Lemma 12. For a sufficiently small e > 0, there exist uniform constants C > 0, x > 0; an d C > 
such that, for every uj, every G (0, 1), and every g G "H_(cr n+m u;; 



(n) 



g o 7* n+m)n+1 di/ n) 



< B_(o~ 



n+m 



m-x log|w| n ^| 



(17) 



/or a// a G -4 n (w) and a// n,m > 0. Here B„(a n+m u;-&-, g) = C||#||oo + K_(a n+m u; 0_, #) tfl 1 
and #(??_) =max(C,/ /2 ) G (0,1). 

We remark that the above lemma is a statement of statistical memory loss starting from an 
initial measure supported on a single unstable curve. This is the strong version of Theorem 8 
hinted at in the last paragraph of Section 4.4. 

Proof of Lemma 12. This follows directly from the time-dependent coupling argument in [26] as 
we now explain. In [26] we proved the analogous statement for "ordinary" Holder continuous 
observables g. Here we just need to check that the estimate one obtains for g G 'H.-{p n+m uj\ is 
the one appearing in (17). Because 



IV 



(n) 



+m,n+l 



dz/ n) 



IV 



(«) 



goJ=, 



n+m,n+l u -'« 



dz/ n) 



9 J 7 ra+m,n+l d/i, 



M 



the desired bound is related to the convergence of the push- forward (J- n +m,n+i)* l/ ct to the (invari- 
ant) measure \i with increasing m in a weak sense. It is shown in [26] that there exists a uniform 
constant x > such that, for m' = [x|log |W+^I|] i the measure (J 7 n +m',n+i)* l 'a is "proper" in 
the sense of [26]. (The terminology originated in [5] and the references therein.) Moreover, the 
measure \i is readily proper. It is therefore certainly enough to show that 



/ 9 ° J~ n+m,n+rn'+l d/i 1 - J go 7, 



n+m,n-\-m'+l 



d/i z 



< Bja 



n+m 



;is) 



for any two proper measures /i 1 and /i 2 . 

We recall from [26] that, for each m > m' + 1, two proper measures // (i = 1, 2) can be decom- 
posed as A*| 1 ( m _ m /)/ 2 ] + Y^j=o~ m ^ flj- Here ti\r m -m')/2~\ are the measures that remain uncoupled 
after \(m — m')/2] steps (counting from time n + ml). Provided e > is small enough, we can 
assume that they satisfy /j, r( m _ m /w 2 ] = A i2 ( m _ m /)/ 2 ] {-M) < CC™ 1 '™" for uniform constants C > 
and ( G (0, 1). Therefore, 



9 -F n+m,n+m'+l d// 1 

< c||dU m -" 



9 3~ n+m,n+m'+l d/i 
f(m-m')/2l 



+ 



J=0 



J 9 J~ n+m,n+m'+l d/ij — J g O J 7 n + m ,n+m'+l d/ij 
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On the other hand, the measures (J r n+m / + j jri+m / + i) ! , t /i* (i = 1,2) are coupled: Both measures 
are supported on a family (called a "magnet") of homogeneous local stable manifolds W^ +m , + j, 
corresponding to the future (u n+m ' + j , u n+m ' + j + i . . . ), in such a way that their masses on each such 
local stable manifold agree. In particular ft^(Ai) = /i|(A^). Notice that, for each pair of points 
x,y in the same local stable manifold W^ +rn , + p 

\g ° J~ n+m,n+m'+j+l(x) ~~ 9 ° J~ n + m ,n+m'+j+l (ij) \ 

< K_{o- n+m u-d_,g)dT m ~ 3 (X ' V) < K_{a n+m u;$„,g)$"- j , 

because g G K-(ff n+m w;t?_). Since ^jftji-M) = 1, it then follows that 



\{m-m')/2\ 

£ 

3=0 



g ° J'n+m^n+m'+l d/ij ~ 9 ° ^~n+m,n+m'+l d/ij 



\{m-m')/2\ 
3=0 

We refer the interested reader to [26] for a detailed construction of the coupling used above. 
Combining the obtained estimates yields (18), which was to be shown. □ 

In order for Lemma 12 to yield a useful bound on the inner integral in (16), we need to assume 
m ^> x|l°g in the exponent appearing in (17). Thus, for a fixed value of m, the partition 

element W+ a should not be too short. The following lemma provides a tail estimate on the 
prevalence of short partition elements. 

Lemma 13. There exists a uniform constant C p > such that the following holds. If u G f2 sa t; 
then the factor measure \W of the measure jj relative to the partition (a;) satisfies 

\ {n) {a G A { ? : \Wi n l\ <£}< C p £ (19) 

for all £ > and for all n > 0. 

It is possible to prove that, for small e, (19) holds for all n G Z and all u> G Q 2 ^sat (starting 
from a generalization of Theorem 5.17 of [6]). However, since the proof is long and technical, we 
resort to the weaker assertion given in Lemma 13, which is perfectly sufficient for our needs. 

Proof of Lemma 13. Denoting Z^ = J^ n {W^l^ 1 d\( n \a) for all n G Z, it can be shown using the 

uniform growth lemma in the time-dependent setting (see [26]) that there exist uniform constants 
C p > and tf p G (0, 1) such that Z^ < |C P (1 + ti$Z®) holds for all i G Z, all n > 0, and 
all sequences uj G VL. Under the saturating past condition (a; G O sat ), u>i = Ui holds for all 
i < i , for some zo, and it is an exercise in the theory of billiards [6] that Z^ ' < oo (because the 
relevant partition corresponds to the homogeneous unstable manifolds of the fixed map Fu iQ )- 
Observing that vf = V { + o) for all i < i , we thus obtain Z^ io+m ^ < |C P (1 + tf^+^Z^) for all 
i < io and all m > 0. Taking the limit i — > — oo, we see in particular that Z^ < ^C p for all n > 0. 
Markov's inequality now yields the result: \^{a G A { + ] : \W^l\ < £} < Z^£ < C p £. □ 
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We are in position to finish the bound (16) on I 2 . For any £ G (0, 1], the decomposition 

4*) = { a G Af : \W^%\ > £} U {a G ^ : \wQ\ < I}, 
in conjunction with Lemmas 12 and 13, yields 

\I 2 \ < sup \f\ ■ B_(a n+m u;$„,g) (0(tf_))™-x|i°g^l . x W{ a e . | W W | > 
+ sup |/| • sup |S| • \ {n) {a G : < £} 

< sup\f\- B_(a n+m u;$„,g) (^(^_)) m ^ x|log ^ + sup |/| • sup \g\ ■ C p £ 

For any m > 0, we set £ = £{m) = e~ m//2x , so that x\ log €| = m/2. Since sup |g| < 2 sup |g|, we get 
the final bound 

\h\ < sup |/| ■ B.^^w; 0_, 0) (^_)) m/2 + 2C P sup |/| sup |<?| e~ m / 2 *. (20) 
Collecting (15) and (20), we arrive at 

f ■ go T n+m d/i - J fdfi J gdfi 

< ll/IU^K^; <?) {6{§_)) m l 2 + 2C P ||/|UMU e~ m ^ + 2 IMU-M"; /) 

Considering separately the two cases m = n and m = n + 1, and replacing n + m by n in the end, 
the bound claimed by the theorem follows. This finishes the proof of Theorem 11. □ 

5.3. Multiple correlation bounds. As mentioned already, the function spaces introduced above 
have a special structure. The following lemmas reveal an important facet of that structure. It 
will serve as the stepping stone from pair correlation bounds to multiple correlation bounds, and 
eventually to the invariance principle. 

Lemma 14. Let tu G £1, $ G (0, 1), and n > 1 be fixed. For all f G "H+(w; $) and g G %-{a n u] ■&), 
we have 

f o jF~ 1 {uj, ■ ) G H + (cr n uj; i?) and g o F n (u, ■ ) G U-{u\ ■&) 

with 

K + (a n u; 0, / o JF-V, ■ )) < 0, /) d n and K„ (w; 0, o JF n (u;, • )) < (<r n u;; 0, g) $ n . 

Proof. Suppose x,y G for some a G -4_. Then J z n (uj,x),J r n (uj,y) G W 7 !"^ for some /3 G A_ , 

which in particular means that s^(a n u; T^ou, x),J- n (u, y)) = n + s_(w; x, y). In other words, for 
geU-(a n u;i9), 

\g(T n (u, x)) - g(T n (oj, y))\ < K_ (a n u; 0, g) $ n d s ~^ . 

Since x, y and a were arbitrary, the claim for g o JF n (oj, ■ ) follows. 

The proof for / o J r ~ 1 (o;, •) is very similar. Suppose x, y G for some a G A { +\ Then 

J 7 " 1 ^, x), J r ~ 1 (w, y) G W+,/3 for some f3 G .4.+ , which implies s + (u; J r ~ 1 (a;, x), J r ~ 1 (w, y)) = n + 
s + ((T n w; x, y). Thus, for / G %+(a>;i?), 

l/^-^a;^))-/^-^^^))! < J^w;*,/)^^^™). 

The lemma now follows. □ 

As a corollary, we get the following result. 
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Lemma 15. Let uj G Q, n > 0, and $o,...$ n G (0,1) be fixed. We set •& = max <j<„ If 
fi G H+^uj; and ||/i||oo < oo /or < i < n, then the function 



/=n^°^ + i(^-) 



i=0 



belongs to % + {a n uj]'d) with 



K + (a n tu;$J) < II H/ilU ) K+iotuirfuMffT. 
i=o / 

Similarly, if gi G 'H-(cr 4 u;; -#i) and ||ft||oo < oo /or < z < t/ien i/je function 

n 
j=0 

belongs to T-L-{u)]'d) with 

i=o Vjyj / 

Proof. Consider some a G .4.+^ and two points x,y £ W+i- Notice that /(x) — /(y) equals 

n / i— 1 \ / n \ 

E nM^I+i^)) [/*(^ifi("^)) -/<(^in(^v))] II bfcj+Mv)) ■ 

i=0 \j=0 J \j=i+l J 

Because ^"^(u;, • ) = Jv^fVw, • ), Lemma 14 yields /j o J 7 "^^, • ) G H + (a n u; with 

ft, /, o • )) < Jjf+^'w; ft, fi) #r- 

Therefore, 



|/(s)-/(y)| < E nn/iii» ) K+i'WuttVr 4 * 

i=0 \j^i 



i 



The proof of the other part is similar, and we omit it. □ 

We are now ready to turn to bounding multiple correlation functions. The following theorem 
gives a rather general and explicit bound. We record it for completeness, as it is interesting in its 
own right. A special case (Corollary 17) admitting a much tidier expression for the upper bound 
will be sufficient for the purpose of proving the main result of the paper. 

Theorem 16. For a sufficiently small e > 0, there exist uniform constants C p > 0, C > 0, £ > ; 
and x > such that the following holds. Let uj G f2 sat , n > 0, m > 0, k > 0, ^ +j0 . . . , G (0, 1) 
and fi^fi . . . , $_ m £ (0, 1) all be fixed. We set $ + = max <i<„^ +i! and = max <j< m As- 
sume that fi G T-L + (a l u; and ||/j||oc < oo for < i < n, as well as g n +k+i G H-(o- n+k+l uj] ft-,i) 
and ||g n +fc+i||oo < oo for < i < m. Denote 

n m 

F = Y[fi° Fi(u, ■) and G = JJ# n+fc+ j o T n+k+i (uj, • ). (21) 

1=0 8=0 
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Then 



FG dfi - / Fd/i / Gdfi 



< 



cH\\g n+ k + i\\oo + E II H^fc+ilU ) K_(a n+k+i u;^, i ,g n+k+i )^i ii ) C 1 

i=0 \i=0 



fc/4 ' " x ' "' 



II ll/<lloo (max(C,/ /2 )) + 2C P H ||/*I|oo ]J \\9n +k+ i\\oo e~ k ^ 

J=0 / \i=0 J \i=0 

m 



+ 2 1 e n ii/iiioo ^+(^5 * /*) n ii^+^iu ) i? 

.i=0 Vj'^j / / \*=o 



fc/2-1 
i II oo | " + 



Proof. Notice the simple identity 

•^(^•^n V, •)) = ^+i( w > •)> < i < n. 
Let us then define the functions 

n 
i=0 

and 

m 

g = G O T^^OJ, • ) = g n +k+i ° ^ n+fc+i.n+fc+ltW; 1 )• 

i=Q 

By invariance of the measure 

J FGdfi = J F{7~ l {uj,x))G{7- l {u,x))d^{x) 

Fi^i", x)) G(J^. k (u, F n+kjn+1 (u, x))) dfi(x) 

f{x)g{J : nJrkyn+1 {u,x))d^{x). 



For the same reason, 

J Fd/i = J /d/i and J G d/i = ^ gd/x, 

so that 

y FGd/i - y Fd/i y Gd/j = J f ■ go JF n+k>n+1 (u;, • ) d/i - J f dfi J gdfx. 

Observing that J- n +k,n+i(Wi ') = J r k(o' n ^ ! ] •)> we wish to apply Theorem 11 next, and thus 
bound the right side of the above expression. To do so, we need to prove that / G H + (cr n w;tf + ) 
and g G "H_(o" fc (<7 n w); $_). By assumption, f\ G 'H + (a l co; for < i < n, and the first claim 
follows immediately from Lemma 15 with 

K + (a n u;# + J) < (nil/ill-) K^juj^Ji)^. (22) 
i=0 ) 
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On the other hand, g = Y\.T=o 9n+k+i ° 

By the assumption of the theorem, g n+k+i G 
T-L-(a' l ((j n+k u)\ for < i < m. Thus, according to Lemma 15, the function g is indeed in 
HJa k (a n uj):$) with 



K_(a k (a n u);^,g) <J2[U Ik+M-ilU ) K^a n+k+i uj-$-^g n+k+l )W 
Now, Theorem 11 yields the bound 



FGdfi- / Fdfi / Gd/i 



(23) 



< 



< 



/ • g o T h {o n u\ -)dfi- J fdfi J gdfi 
B„(a n+k cu; 0_,g) {d{§^) k " + 2C P ||/|| 00 ||</|| 00 e"*^ + 2K + (a n co; + , /) |MU tff" 1 

C7lJ||^ rH . fc+i || 00 + if_( < 7 n+fc a;;i?_^)C 1 [ (0(tf-)) fc/4 



nn/*ii' 

,i=0 

+ 2C p (ni|/ 1 



8=0 



v i=0 



n+fc+i II oo 



n+fc+i II oo 



I? 



fc/2-1 



v« = 



1 /2 

Inserting (22), (23), and 0(i?_) = max(C,0Z ) G (0,1) proves the theorem 



□ 



Now, we specialize to the actual case of interest, which is related to proving convergence of the 
characteristic function of the vector- valued process in Theorem 1. Thus, below, the functions f\ 
and g n +k+i are going to be uniformly bounded by 1, as well as have the same rate and the same 
constant in the dynamical Holder continuity conditions, as follows. 

Corollary 17. For a sufficiently small e > 0, there exist uniform constants C p > 0, C > 0, ( > 0, 
and x > such that the following holds. Let u> G fi S at; n > 0, m > 0, k > 0, $ G (0, 1), and K > 
be fixed. Assume that fa G H + (o- l u;'d) with K + (o- l u;$, fi) < K and ||/j||oo < 1 f or < i < n, 
as well as g n+k+i G T-i^(a n+k+t u] ■&) with K_(a n+k+l u; g n+k+i ) < K and \\g n+k +i\\oo < 1 for 
< i < m. With F and G as in (21), we have the bound 



|Cov„ (F,G)\ 



< AX 



FGdfx- J Fdfx J Gd/i 

where A = C + 2C P + 3K(1 - tf)- 1 ^- 1 > and X = (max(C, 1/2 , e" 1 ^)) 1/4 G (0, 1) . 
Proof. From Theorem 16, we easily obtain 

( n "\ n 

|Cov M (F,G)| < Ic + K^&ti- 1 I (max(C,^ 1/2 )) fc/4 + 2C p e~ fc/4x + 2ir^?9 n -^ fc/2 - 1 . 



i=0 



i=0 



The result follows from this, because ^27=0^ = SILo^™ * — (1 — ^) 



□ 



Notice that the above multiple correlation bounds were deduced directly from the pair correlation 
bound in Theorem 11. 

5.4. Characteristic function bounds. Let us emphasize that our pair and multiple correlation 
bounds hold for all sequences u G f2 sa t, without any statistics required on the sequences. We will 
show below that an application of the bounds on mixing, asymptotically mean stationary, random 
sequences yields the vector- valued almost sure invariance principle. 

The next theorem follows when [8] is applied to the special case of bounded processes. 
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Theorem 18 (Gouezel [8]). Let d be a positive integer and consider a sequence (A n ) n > ofR d - 
valued random variables which is centered and bounded. Given integers n > 0, m > 0, < b\ < 
62 < • • • < b n+m+ i, k > 0, and vectors ti, . . . , t n+m G M. d , set 

m bj+i—l+k 

*£i = £v E A * ( 24 ) 

j=n i=bj+k 



C{n+m) 

(25) 



for brevity. Now, suppose there exist constants t > 0, C > 0, and c > such that 
E (V<n+i*Si,„ +m ) _ E (^x^VjE {e lX -+^+^j < Ce- Ch M + ^max \b j+1 - b 

holds for all choices of the numbers n, m, bj, k > 0, and all vectors tj satisfying \tj\ < t. 

(I) If (A n ) n > is stationary, there exists such a symmetric, semi-positive-definite, dx d matrix S 2 
that the statements (2)-(4) listed as consequences in Theorem 1 hold true. 

(II) If (A n ) n >o is non-stationary, and if there exists such a d x d matrix S 2 that, for any a > 0, 



sup n 

£>0,n>l 



< 00 (26) 



te+n-i e+n-i \ 
E A fe ® A fe - n5]2 
k=i k=£ / 

is satisfied, then items (3) and (4) in Theorem 1 hold true. (Here || • || denotes matrix norm.) 

Denote Q + = (B E (0)) N and let P be a probability distribution on it, as in the statement of 
Theorem 1. Notice that we can embed fl + in f2 sat canonically: Given a sequence u G fi+, we 
identify it with the sequence Co G Q SSLt given by u)j = uji for all i > and u)j = u for all z < 0. 
For this reason, we will not make an explicit distinction with one- and two-sided sequences in 
the following argument. (That is, u> G Q + is picked randomly according to P and fixed, after 
which it is augmented to an element of f2 sat , still denoted co, which is then used in the subsequent 
computations.) 

In our application of Theorem 18, given oj G Q + , we set A n (uj,x) = f (u n , J n (w, x)) and define 
Xn)n(uj,x) as in (24). Notice that we can write 

b n +l— 1 

F(u,x) = exp(iX[°l(u,x)j = JJ fi(u i} ^{u, x)), 

i=b\ 

where fi(u>i, x) = exp(itj • f (cji, x)) with j G {1, . . . , n} chosen so that bj < i < bj + ±. Similarly, 
G(u,x) = exp(iX ( n % n+m (uj,x)) = \ \ g^Ui, ^(u, x)), 



i=b„ +1 +k 



where g^u^x) = exp(itj • f(ui,x)) with j G {n+1, . . . ,n + m} chosen so that bj + k <i< bj + i + k. 
Thus, our goal is to show that the covariance of F and G with respect to the measure P <g) \i has 
a uniform upper bound like the right side of (25). Notice already that 



Cov P0M (F, G) = j Cov^ (F, G) dP + Cov P (J F d/i , J G d^ . 



(27) 



Both (uj,x) 1 — y f(u n ,x) and (u,x) >->■ F n {u,x) are measurable for every n > 0. (The first one 
of the two is an assumption, and the second is a consequence of the continuity statement after 
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Lemma 5.) In particular, (uj,x) (-)■ F(u,x) and (uj,x) !->■ G(u,x) are measurable and absolutely 
bounded by 1. Thus, all the integrals (including the covariances) in (27) are defined. 

Owing to the Holder continuity assumption on f , Lemma 9 shows that there exist constants 
i? G (0, 1) and K > such that, for every i > 0, every vector component of the function f (cUi, ■ ) 
belongs to H-(a l u; fl V. + (a t u; and 

\i{ui,x) - f(ui,y)\ < M s±{cTl ^ y) Vx,ye W^ a Va G A ( £ . 

(In each condition we either choose "+" everywhere or "— " everywhere.) Assuming |tj| < t for all 
j > 1, the bound \e m — e lb \ < \a — b\ for a, b G IR implies 

\fi(ui,x) - fi(u i} y)\, \gi(ui,x) - 9i{^y)\ < t |f(w<,a;) - f(u>i,y)\. 

We can thus apply the multiple correlation bound of Corollary 17, defining K = tK. The small 
formal difference that here the product in the expression of F starts from time b\ instead of is 
abolished by regarding momentarily /j = 1 for each i = 0, . . . , b%— 1. Consequently, | Cov M (F, G) \ < 
AX k , with the uniform, non-random, constants A and A. In particular, 



Cov^ (F, G) dP 



< AX k . 



The maps F : u !->■ J -F(w, • ) d/i and (5 : u M- / G(a;, ■ ) d/x are measurable. What is more, 
F{oj) only depends on Ui with bi < i < b n+ i, while (5 (a;) only depends on Ui with + < 
i < b n+m+ i — 1 + k. Hence, if the sequence (cjj)j> is random and with < a < b integers, 
stands for the sigma-algebra generated by {ui)\ =a , then F(oS) and G{oS) are, respectively, $ n+1 1 - 
and 5^ +i+fc -measurable complex- valued random variables. By the rho-mixing assumption (Al), 
and the fact that \F\ < 1 and \G\ < 1, the identity in (5) yields 

|Cov P (F,C7)| < ^(jfe) < Be~ ck . 

By (27), we conclude that, in the notation of Theorem 18, 

|Cov P ^ (F, G)\ < AX k + Be~ ck 

holds uniformly for all choices of the numbers n, m, bj, k > 0, and of the vectors tj satisfying 
|tj| < t. The obtained bound is obviously stronger than Gouezel's condition (25). Hence, the proof 
of Theorem 1 in the stationary case is complete. 



5.5. Convergence to limiting covariance. In this section we establish uniform convergence to 
a limiting covariance matrix, which is needed in the non-stationary case; see (26). We denote 

e+n-i 

S n ,i = Afc and S„ = S n) o- 

k=l 

Lemma 19. The matrix S 2 in (9) is well defined. It satisfies the bound in (26), for any a > 0. 
Proof. By /x-invariance, 



E(S n ,, <g> S n , e ) = / / (S„ ® S n ) (cr e w, x) dfx(x) d¥(u). 
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/ (A fc ® A fc ) (w, • ) d// + / (Aj ®A k + A k ® Aj) (w, • ) d/i 

fc=0 0<j<fc<n-l ^ 

n— 1 „ „ 

^ / (A ® A ) (a fc u;, ■ ) d/i + / (A <g> A fc _,- + A^ <g> A ) (<A;, • ) d/i 



fe=0 
n-l 



0<j<fc<n-l 
n— 1 n— 1— m 



Denoting 



and 



lb J_ n r (■ J_/t _L 1 1 L n 

/ (A <8> A ) (^'w, ■ ) d/i + / (A ® A m + A m <g> A ) (<r J 'u;, • ) d/i. 

j=0 m=l j=0 

V (w) = y"(A ® Ao)(w,x) d/i(x) 
V m (w) = y (A <g> A m + A m ® Aq) (w, z) d//( 



x), m > 1, 



m=0 



n-l 



as well as using the notation introduced in (6), we thus have 

n— 1 n— 1— m n— 1 

E(S rtj/ (8) S n/ ) = ^ E(V m o^) = ^(n-m)(V ra o4_, 

m=0 j=0 
n-l 

n 

m=0 m=0 

By Theorem 11, there exist C > and /3 > such that 

sup ||V m (w)|| < Ce-^ m , m > 0. 

Therefore, ||(V m o <r £ )fc|| < Ce _/3m for all m,£ >0 and > 1, which yields 



^ (V m o cr £ ) n _ m - m (V m o o- £ ), 



(2f 



n-l 



m (V m o a e ) n „ 



m=0 



n-l 



m=0 



m=0 



or 



n-l 



E(S n/ ® S v ) - n (V m o </) r 



m=0 



< C", I > 0, n > 1. 



(29) 



Notice that V m is bounded and S^-measurable. By assumption (A2), there exist the limits 

(V m )oo = lim (V m o a l ) k , I > 0, m > 0. 

k— >OD 

Recall that in (9) we have claimed that the desired covariance matrix is given by 

oo 

S 2 EE ( V ™)oo- (30) 

m=0 

Certainly by (28) the series on the right converges in norm and yields a well-defined matrix. 
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Lemma 20. For any bounded and 5 -measurable function g : (i? e (0)) N — > R d ; 

(g°A» = (g>oo, ^>0, m>0. 

Proof. We only check the claim for t = 1, the general case being similar. Clearly g o a is bounded 
and measurable (5 being the product Borel sigma-algebra) , and 

(g o a) k - (g) k = - (E(g o a) - E(g)) ^ as k -»■ oo. 
Since (g o er) fe -y (go er)^, the claim follows. □ 
Continuing with (29), if \ogn < n, 

n— 1 /3 _1 logn— 1 ra— 1 



m=0 



m=0 



m=/3 — ilog n 



where 



n-l 



n-l 



m=/3 _1 log n 



E (V m o 

m=/3 _1 log n 

In other words, we have shown that 

/3 logri— 1 

E(S n/ g> S n ,*) - n ^ ( V ™ ° a& )n-m 



-/3m < ^ ( 

~~ n 

m=0 



-j3m 



m=0 



<2C", f3~ l \ogn<n. 



Next, recalling (30) and Lemma 20, we decompose 

/3 _1 logn— 1 oo /3 — logn— 1 

£ (V m o4_ m = S 2 - >T (V m ) 00 + £ ((V m oa £ ) n _ m -(V m ) c 

m=0 m=f}- 1 \ogn m=0 

The middle term on the right can be bounded using again (28). Indeed, 

C 



m=/3~ 1 log n 



m=/3 1 logri 



n 



Hence, using assumption (A2), we can compute 

/8- 1 log n-l 

|E(S„,* ® S„, € ) - nS 2 || < 3C7' + n ' 



E ° Cr e )n-m - (V m ) c 



m=0 



/3 ^ogn— 1 



/3 Uogn— 1 



< 3C" + n E G ' 



/» ' n—m » m oo 



< 



3C" + C7n E C < 



-/3m 



m ' n—m ' 



m=0 



m=0 



< 3C" + OD^ log n, /3~ l \ogn<n. 
Finally, (29) yields also the crude bound 



|E(S n/ <g> S n/ ) - nS 2 || < C + n 



n-l 



E ( V m o CT e )n-rr 



m=0 



< C" + n(||S 2 || + C) 



< 2C'(1 + /rMogn), /TMogn^n. 
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Collecting the bounds in the two regimes (3 1 logn < n and /3 1 \ogn > n, we see that there exists 
a constant C" > such that 

||E(S n g> S n ) - n£ 2 || < 3C" + C'logn, n > 1. 

This finishes the proof of Lemma 19. □ 

An application of Lemma 19 with part (II) of Theorem 18 finishes the proof of Theorem 1. □ 

Remark 21. Related to the special case of Theorem 4, the above proof contains the interesting 
fact that a pair correlation bound (for a well-chosen class of dynamically Holder observables) alone 
implies the vector-valued almost sure invariance principle for Sinai billiards with fixed scatterers. 
This extends the analogous result of [24] about the central limit theorem for Sinai billiards. Before 
that, the otherwise classical central limit theorem for Sinai billiards had been obtained via the new 
method of multiple (as opposed to just pair) correlation bounds in [5, 6]. The author has learned that 
also Pene [19] has proved the central limit theorem for Sinai billiards using correlation functions 
without, however, pursuing finer limit theorems. 

5.6. Proof of Lemma 2. Let us write F k for the probability measure on ((i? e (0)) N , 5) given 
by F k {A) = f Ej; lp O~^) for a11 A e & % assumption (A2), F k {A) = (l A ) k (l^ as 
k — > oo. The Vitali-Hahn-Saks theorem now states that the map P : $ — > [0, 1) : A (l^)oo is a 
probability measure. The measure P is called the stationary mean of P, which is justified by the 
simple but crucial observation made in Lemma 20 of the preceding section that P — unlike P - 
is invariant for the left shift a. Because /i is invariant for each of the billiard maps, the measure 
p = Pcg)^ is therefore invariant for the skew-product map $ defined in the statement of the lemma. 
Observe that, writing E for the expectation relative to P, (9) reads 

oo 

S 2 = E(A (8) A ) + B( A o ® A ™ + A m <g> A ). (31) 

m=l 

Here A n (u,x) = f(u n ,J- n (u,x)) = f ($ n (o;, x)) with the understanding that f(u,x) = f(uj ,x). 
The following is a complementary result to Lemma 19. It compares the covariance of the sum 
S n = Y17=o A * = J2*i=o f ° < ^ )i ' computed with respect to P instead of P, to nS 2 . Taking advantage 
of the $-invariance of P yields a match better than what one would naively infer from (26): 

Lemma 22. There exists a constant C > such that 

sup||E(S„(g)S n ) -«S 2 || < C. (32) 

n>0 

Proof. We start with the elementary identity 

n— 1 

E(S n ® S n ) = n E(A ® A ) + (n - m) E(A ® A m + A m <g> A ), 

m=l 

which is where we need invariance. Recalling (31), we have 

oo 

E(S n ® S n ) - nS 2 = a n( m ) E(A ® A m + A m ® A ) 

m=l 

where a n (m) = —m for 1 < m < n and a n (m) = —n for m > n. Because |a n (m)| < m, 

oo 

||E(S„® S„) -nS 2 || < ^m||E(A ® A m + A m ® A )||. 

m=l 
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In the notation of Section 5.5, 



j fc-i 

E(A <g> A m + A m ® A ) = Km - E(V m o a 

fc^oo A; z — ' 



£=0 



The uniform bound in (28) yields ||E(V m o a e )\\ < Ce /3m and hence the result. 



□ 



Recall the /z-average of f vanishes identically. Given a vector v G M. d , we define f v = v T f . Since 
v T (A m <g> A n )v = v T (A„, ® A m )v = / v o $ n ■ / v o $ m for all m, n > 0, (32) gives 



n— 1 



Var P ^ / v o $ fc - n v T £ 2 v 



vfc=0 



|v T (E(S n <8> S n ) - n£ 2 )v| < C"|v| 



(33) 



uniformly for n > 1. 

Suppose E 2 is degenerate. In other words, there exists a vector v G M d such that S 2 v = 0. By 
the bound above, VarpQ2£~ ^ v ° ® k ) — C'l v | 2 uniformly. Owing to the $-invariance of P, we can 
therefore apply [13, Lemma 1] (see also [21]) and conclude that f v must be an L 2 (P)-coboundary. 
In other words, there exists g 6 £ 2 (P) such that f v = g — go& holds P-almost-everywhere. Since 



j=0 J 



to, x)\ 2 dP(a;, x), 



(34) 



the claim of the lemma in one direction follows. 

To prove the claim in the other direction, suppose f v = g — go§ for some g that satisfies the 
condition given in the lemma. Then g G L 2 (P) by (34), and ^"Tq 1 / v o $' = g - g o $™ holds 
P-almost-everywhere. This immediately gives 



'n-l 



Var P o $ fe = \\g - g o $ n \\ 2 L2{P) < M\g\ 



Combining the bound with (33) we get v T S 2 v < n 1 (C|v| 2 + 4||(7|| 2 2( .p^) for all n > 1, which is 
only possible if S 2 is degenerate in the direction of v. This completes the proof of Lemma 2. □ 
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